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Abstract 

The small polaron with generic, nondiagonal boundary terms is investigated within the 



1 

rrt framework of quantum integrability. The eigenvalues of the model are extracted by using the 

I fusion hierarchy of the transfer matrices and the eorrespond.ng Bethe Ansat. equations are 

presented. For particular values of the anisotropy parameter the fusion hierarchy truncates, 

^ giving rise to a set of functional relations for the transfer matrix. Exploiting the latter ones, 

i/~J the same set of eigenvalues is rederived, confirming our results. Finally, we comment on 

CN the eigenvectors of the model and explicitly compute the state with all sites unoccupied for 

^ arbitrary chain lengths. 

m 
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1 Introduction 

The investigation of exactly solvable models over the last many years has been proven very 
fruitful in extracting physical information regarding phenomena lying in strongly coupled 
regimes [1]. For integrable models in one spatial dimension, the framework of quantum 
inverse scattering method (QISM) [2, 3] provides an efficient way of computing the various 
physical quantities. In the present article we study the small polaron model, which can be 
also regarded as a graded version of the XXZ quantum spin chain. The i?-matrix satisfies 
the graded Yang-Baxter equation [4, 5, 6] 

Rl2{\) Rl3{\ + /i) i?23(/i) = i?23(/i) Rl3{\ + ^l) RuW , (LI) 

and the resulting bulk Hamiltonian coincides with the bulk Hamiltonian of the XXZ spin 
chain, after employing a Jordan- Wigner transformation. The effects of supersymmetry be- 
come visible when periodicity of the boundary conditions is relaxed, hence more general 
boundaries are considered, corresponding to open boundaries in the spin chain picture. 

Within a graded version of Sklyanin's reflection algebra [7] these boundary conditions can 
be implemented while keeping the integrability of the model. For purely diagonal boundary 
flelds, i.e. respecting the f/(l) symmetry of the bulk, the spectrum of the model can be 
obtained using algebraic Bethe Ansatz (ABA) methods. For non diagonal boundary condi- 
tions, however, this symmetry is broken and a simple reference state does not exist anymore, 
rendering the ABA framework insufficient. For the respective ungraded model various meth- 
ods have been succesfully employed in the past, such as the fusion hierarchy of the transfer 
matrices and its truncation [9, 10], the construction of a vacuum state by using gauge trans- 
formations [11], generalized TQ equations [12] as well as functional relations derived directly 
from the Yang-Baxter algebra [13] and by separation of variables [14, 15, 16] amongst oth- 
ers. Nevertheless, in all of these approaches the actual solution of the eigenvalue problem 
relies on the boundary parameters satisfying certain algebraic relations or is restricted to 
particular values of the bulk anisotropy parameter. 

Motivated by a recent analysis of the free fermion case [17] where the Grassmannian 
nature of the nondiagonal boundary parameters was proven sufficient to solve the model 
without imposing any constraints, we studied the small polaron model with nondiagonal 
boundary conditions, expecting that a similar situation may also hold. As it will become 
transparent in the next sections, it turns out that this is the case for the interacting fermions 
as well: supersymmetry lifts the need of imposing constraints on the boundary parameters. 
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Furthermore, the structure of the eigenvectors is greatly restricted, so that certain eigenstates 
can be computed exactly for an arbitrary number of chain sites. 

The paper is organized as follows: first, we describe the basic facts of the small polaron 
model and setup our conventions. In Section 3 we study the eigenvalue problem for the 
transfer matrix (and consequently the Hamiltonian) with non-diagonal boundary conditions 
based on the fusion hierarchy obtained in [18]. Assuming the existence of a particular 
limit the spectral problem can be formulated as a functional equation which in turn can 
be solved by deformation of the corresponding problem in the diagonal case. The transfer 
matrix eigenvalues are found to depend on two distinct sets of Bethe roots, which satisfy 
two coupled sets of Bethe Ansatz Equations (BAE). In Section 4 we study the same problem 
restricted to particular values of the anisotropy parameter, for which truncation identities 
of the fusion hierarchy were found [18]. The functional relations for the transfer matrix 
are derived and written into a vanishing determinant representation, which leads to the TQ 
equation without the need of additional assumptions. The eigenvalues of the transfer matrix 
coincide exactly with those obtained directly from the fusion hierarchy, thus ensuring the 
validity of our results. In Section 5 we compute the 'vacuum eigenstate' of the model, i.e. 
the unique state which reduces to the Fock vacuum in the limit of diagonal boundary fields. 
Based on this result we propose an expression for the generic structure of the eigenvectors 
of the model which allows for a complete solution of the spectral problem in principle. We 
conclude with discussing our results and future directions. 

2 The small polaron with open boundary conditions 

The small polaron model [19, 20] describes the motion of an additional electron in a polar 
crystal. The physics of these interactions is captured by the following bulk Hamiltonian 

-f^buik = X^ -^ (cj+1 Cj + c] Cj+i) + V [rij+i rij + Uj+i Uj) , (2.1) 

where c|. and Ck denote the creation and annihilation operators of spinless fermions at site k 
respectively, obeying anticommutation relations {c^,q} = 5ki- We have also introduced the 
occupation number operators n^ = clck and n^ = 1 — Tifc, so that the parameters t and V may 
be interpreted as hopping amplitude and density-density interaction strength respectively. 
It is possible to derive the above Hamiltonian through the QISM framework, thus rendering 
the system integrable [21]. In the Hamiltonian (2.1), periodic boundary conditions are to 



be assumed. However, periodicity can be relaxed and one may consider integrable open 
boundary conditions, by using the framework of quantum integrability for models with open 
boundaries [7]. The full Hamiltonian of the system contains the additional boundary terms, 
which in general may be non-diagonal ones 



H — -f^bulk + -f^diag + Hj 



nondiag 



Since we deal with a fermionic lattice model, the local space of states is Z2-graded [23]. The 
tensor product is graded according to the rule 



(A 0s Bf\^ = (_i)b(«)+p(f')]pW^«^5^^ 



d 5 



(2.2) 



where the parity p{a) is equal to zero (one) for bosonic (fermionic) indices. All matrix 
operations below, such as the super trace of a matrix, are then to be understood as operations 
on super matrices. We omit the definitions of these matrix operations and refer the interested 
reader to the references [22, 23], whose conventions we follow. 

The fundamental super transfer matrix of the model is given by the super trace 

t{u) = stro\ K+{u) T{u) K~{u) f{u) 



(2.3) 
where the monodromy matrices are defined as 

T{u) = Rmiu) ■ ■ ■ R2o{u) Rio{u) , f (iz) = RoK-u) R^^'i-u) ■ ■ ■ Rmi-u) . (2.4) 



The i?-matrix is given by 



R(u) 



/sin(-u + 2ri) 



sm{2ri) 









sin(-u) sin{2ri) 

sm{2ri) sin(-u) 

\ -sm{u + 27])J 



(2.5) 



satisfies the graded Yang-Baxter equation and enjoys several useful properties, such as uni- 
tarity 

Ri2iu)R2ii-u) = C{u) , (2.6) 



crossing symmetry 



and periodicity 



R',{^i-u-Ar])R',\^{u) = au + 21^), 
Ri2{u + tt) = -alRi2{u)al = -alRi2{u)al . 



(2.7) 
(2.8) 



In the above we have also defined 

au)=giu)gi-u), and giu) = ) ,^ ' ■ 2.9 

The i^-matrices, which contain the boundary information, satisfy the graded refiection 
algebra [7, 8] and have the following generic expressions 

/sin(-u + ■?/'_) a" sin(2-u) 

K (u) = LU I 

/3" sin(2-u) — sin (m — ■?/'_) 

(2.10) 

^^,, , ,, sin (-u + 2?7 + ■?/'+) a+ sin(2M + 4?7) 

A (M) = u I 

/3+ sin(2-u + Arj) sin{u + 27] — ■?/'+) 

with normalizations cu* defined by 

cj- = and UJ+ = . (2.11) 

sin-?/'. 2 cos 2?] sin ■?/'_i_ 

The boundary parameters ip^ are even Grassmann numbers with a non-vanishing complex 
part, while the parameters a*,/3* are odd Grassmann numbers, being subject to the con- 
dition a^ ■ P"^ = 0. Through the framework of QISM, the super transfer matrix (2.3) gives 
rise to a commutative family of conserved charges, among them lies the Hamiltonian of the 
small polaron with open boundaries, which can be derived as 

d 



*'<"' 



= 2iJ + const. (2.12) 

u=0 



and coincides with the expressions (2.1), (2), after identifying t = —esc 2?] and V = cot 2?]. 
The boundary terms read as 

-f^diag = Af+flN -J^-riN + - cot ^_ {ni - 111) 

-f^nonddiag = CSC ^_ (ft^Ci - f3' c{) + CSC V^+ (ft+CAT - /3+cJv) , (2.13) 

where the following shorthands have been also introduced: A/± = ^ csc(2r7) csc^ip^) sin(2?7 ± 
^+). 

The super transfer matrix enjoys crossing symmetry 

and periodicity 

t{u + n)=t{u). (2.15) 



For future use, we also note that the open transfer matrix is normahzed as t(0) = 1 and 
becomes diagonal in the semi-classical limit r^ — > 

(_l)N N 

t(u) = — ^—^ — - 2 sin^ u cos^ u (/?+«- - «+/?-) Ff af 
^ \=o smV^+smV^_L ^ II '^ 

— ( cos^-u sin ■?/'_ sin ■?/'+ + sin^M cos ■?/'_ cos ■?/'_i_) 1 . (2.16) 

Finally, taking the limit 2; = e*" — > oo, the asymptotic behavior of the super transfer matrix 
is obtained, which is needed for later comparison. The leading term contains only odd 
Grassmann boundary parameters and has the expression 

+ - ^ 

t{z) = {-If'^^ z^e^''' (/3+a- - «+/?-) JJ(n, - e^^^OK' + e'^^O + ^i^^) ■ (2-17) 

i=i 

This particular combination of the odd Grassmann boundary parameters emerges in many 

different relations, and as will become transparent below it appears uniquely in all eigenvalues 

of the transfer matrix, thus henceforth it will be denoted as 

g = /3+a--a+/3-. (2.18) 

3 Bethe Ansatz solution from the fusion hierarchy 

In order to solve the model under consideration, we exploit the fusion hierarchy of the transfer 
matrices. Having a transfer matrix at hand, one may construct a family of commuting 
transfer matrices derived from auxiliary spaces of higher dimension, through a suitable fusion 
procedure [25, 26]. The super quantum determinant of the small polaron is given by 

A{u) = { .t ^ ) g{-2u-6r])g{2u + 2r])detK+{u)detK-{u + 2r]). (3.1) 

Note that this expression coincides with the quantum determinant of the model with diagonal 
boundaries as a consequence of the constraint a* ■ /3* = for the off-diagonal elements of 
the boundary matrices. It is useful to introduce the rescaled quantities 

Hu) ^ ^(2^ + 47^) ' ^^^ i^''\n) ^ - (jl^r^u)] &\u) , (3.2) 

with t(°) = -t(°) = 1 and t^^) = t^^\u) = -t{u) and 

n 

e„(M) = Yl C(2m +[n + k]- 2r]) . (3.3) 

fc=i 



In terms of these rescaled quantities, the fusion hierarchy reads as [18] 

?"'(m) ■ i^'\u + n-2r]) = i^"+'\u) - A{u + [n - I] ■ 2r]) ■ i^"-'\u) . (3.4) 

Since the fused transfer matrices commute with each other, the fusion hierarchy can be 
also read as a set of relations between their eigenvalues, which after shifting u -^ u — n ■ 2r] 

gives 

~, A("+i)(M-n-2r7) ~ A("-i)(m - n ■ 2r7) , , 

A{u) = ^—^ --^{u-2r])^—^ ^. 3.5 

AH(M-n-2r7) M'^) (u - n ■ 2ri) 

Assuming that the limit Q{u) = lim„^oo A*^"^(-u — n ■ 2r]) exists then [27], one obtains the 
formulation of the spectral problem as a TQ-equation from the fusion hierarchy 

A{u) = Ha{u) — -=^— Hs{u) ~ — , (3.6) 

Q[u) Q[u) 

the sole requirement being that the //-functions should factorize the quantum determinant 
as Hs{u - 2r])Ha{u) = C-\2u)A{u - 2r]). 

Since the diagonal case should be contained within the construction as a special limit, 
we propose the following deformations of the iZ-functions 

Ha,5{u) = K^s{u) {l + Q fa,5{u)) , (3.7) 

where ha,5{u) are the /i-functions of the diagonal case, which essentially factorize the quantum 
determinant 

, / \ + _ sin(2-u + 4?7) . , , . . , , . f - sm(u + 2r])\ 

ha{u) = u^u —— — — — sm(M + V^+) sm(M + ^_) —— -— 

sm(2-u + 21]) \ sm(-u — 2?]) J 

sm 2u , . . , . ^ / — sm u ^ 



^^(") = ^'^\in(2« + 2r,) ^^"(" + ^r^ " ^^^ ^^"(^- " ^ " ^r^^ ^^.^^^ ^ ^^^ ^.^^^ _ ^^^ 

(3.8) 
and fa,5{u) are generic functions to be determined. There exist two constraints on the 
expressions of these unknown functions. First, their contribution in the factorization of the 
quantum determinant should vanish, yielding 

fs{u-2r]) = -fM- (3.9) 

Furthermore, assuming that the Q-functions have crossing symmetry, the if -functions are 
related through (^^(-u) if Q,(-u) = —(^{—u—2ri)Hs{—u—2ri). This relation holds automatically 
for the /i-functions of the diagonal case, whereas it also provides a second constraint on 
fa,s{u): 

fsi-u-2r]) = Uu)- (3.10) 
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A set of solutions to these functional relations is given by 

fg(u) = Wsm{2u + 4:r]) 

f^(u) = ->Vsin(2M), (3.11) 

with the coefficient W to be determined from the asymptotic behavior (2.17) of the transfer 
matrix eigenvalues. It should be stressed out that the above choice is not unique, but it was 
found to be the only one consistent with the various limits and with the constraints arising 
from the functional relations, as they will be described in the following section. 

For the nondiagonal case that we investigate, a similar deformation should be considered 
for the Q-functions as well, since they are assumed to originate from the higher spin eigen- 
values through a limiting procedure. The most general ansatz extending the Q-functions of 
the diagonal case, Q{u), and containing the special combination of Grassmann numbers, Q, 
reads as 

M 

Q{u) = Q{u) + g b{u) , Q = Y[sm{u-vf'^)sm{u + vf^ + 2r]), (3.12) 

e=i 

with b{u) being a function to be determined. Substituting the deformations (3.7) and (3.12) 

into the TQ equation, the latter one becomes 



A(.)^A--(.)(l-^|i)+^ 






.,Hu+^ Q{u + 27]) 

-hs{u) h fs{u) 



Q{u) ' "''^' Q{u) J\ ' ^^-^^^ 

with the eigenvalues of the diagonal case being given by the relation 

A-'-(") = '=«(")%^ - l^siu)'^^^. (3.14) 

Q[u) Q[u) 

Note that only the Q-functions of the diagonal case appear in the denominator of A(-u). 
Therefore, their zeros v\ are determined by requiring analyticity of the eigenvalues of the 
transfer matrix, in particular of their diagonal part A'^^'^^{u). This yields the BAE 

/i«(t;f) Q(t;f + 2r/) 



,(0)^ nf.So) 



hsivn QK'-^v) 



(3.15) 



or upon substitution 



sin(f) + 21]) \ sm{Vf + ■?/'_) sin(f^ + ip+) 



sint>) J sin(t>) + 21] — ip^) sin(t>) + 2i] — ^_) 



-p-j- sin(f) Vj + 2i]) sin(f£ + Vj + Aij) 
j-J sin(t;f-t;f -2r^)sin(^;f +^;f) 



In the case of nondiagonal boundaries, additional requirements are needed to determine 
the function b{u) in (3.12). Again the choice of this function has to guarantee the analyticity 
of the additional terms appearing in the functional equation (3.13) for the eigenvalues. In ad- 
dition, b{u) should enjoy crossing symmetry. In this spirit, we propose that the nondiagonal 
correction to the Q-functions is given by 

M' 

b{u) = Yl sin(M - t;f ^) sin(M + v^l^ + 2ri) , (3.17) 

with {v^ } 7^ {f) } in general. Matching of the asymptotic behavior dictates that the upper 
limits of the products in eqs. (3.12) and (3.17) should be equal, M' = M. Furthermore, the 
asymptotics (2.17) provide the explicit expression for the coefficient W as 

" sin[^+ + ij_ + {N-2M- l)2r]] ' ^^'^^^ 

Apart from the asymptotics, one should also consider additional limits of the TQ equation 
in order to check the consistency of the procedure. The limit -u — )■ gives A(0) = 1, which 
stems from the normalization of the transfer matrix. Setting -u = in the TQ equation, after 
a quick inspection one is lead to the following constraint for b{u) 

b{-2ri) = b{0) , (3.19) 

which is a special case of a crossing symmetry requirement and automatically satisfied by 
the choice (3.17). Finally, it is interesting to consider the semi-classical limit r^ — )■ 0. In this 
limit the b{u) functions drop out from the TQ equation and the eigenvalues become 

A{u) = {hM - hs{u)) + g fM {hM + hs{u)) , (3.20) 

finding perfect agreement with the semi-classical limit obtained from the transfer matrix 
(2.16), after recalling the definitions (3.8), the solutions (3.11) and the expression (3.18) of 
the coefficient W for r^ — )■ 0. 

The parameters v^ in (3.17) are determined by requiring analyticity of the eigenvalues 
of the transfer matrix for the nondiagonal case. The purely complex part of the TQ equation 
(3.13) yields the already derived BAE (3.16), while the one containing Grassmann numbers 
provides a second set of relations, involving both sets of roots v^^\ v^^\ which reads as 

haivf^^ Adi'^g(t;f)6(t;f) &(t;f + 27^) + ^(^f )g(t;f + 277) 

hsivf) h,{vf) {b{vf^ - 2r^) + Uvf^)Q{vf - 2r^)) ^ b{vf^ - 2r^) + Uvf)Q{vf - 2^) 

(3.21) 



Note that the roots v'^^^ and v^^^ also enter this set of equations through the functions Q{u) 
and b{u). 

In summary, we have two sets of algebraic equations fixing the parameters appearing 
in the ansatz for Q{u). These equations are of nested Bethe Ansatz type, similar to those 
appearing in integrable models based on higher rank symmetries. 

In order to further test our proposed expressions for the eigenvalues as computed from 
the TQ equation, we have analytically diagonalized the transfer matrix for a small number of 
chain lengths. The comparison of the exact eigenvalues with the ones obtained from (3.13) 
exhibits perfect agreement. 

4 Truncation of the fusion hierarchy and functional re- 
lations 

The findings of the previous section are based on the assumption that the limit Q{u) = 
lim„_>.oo A*^") (m — n ■ 2r]) exists. In order to enforce the validity of our results, we follow 
here a different path. For particular values of the quasi-classical parameter, r]n = ^r^, the 
fusion hierarchy truncates and the derived functional relations can be used to derive the TQ 
equation. The truncation identites have been worked out in [18], where it was found that^ 

t»(w,r^„_i) = 0jf_i(w) ■ I - <l>l_,{u) ■ i^-~^\u + 27^„_i,7^„_i) , (4.1) 

where the functions 0jf (m), 0^('u) are given by the following expressions 

<Pniu) = Ml^'iu) /i:(M)/i„(M)[z/:(-M)z/„(M) + U^{u)p^{~u)] 

cl>l{u)=C''{u)^l{u)^-{u), (4.2) 



with 

Mniu) 



/^n(w) 



v^iu) 



1/2 \"sin([n+l]M) 



sin(2?7„) T-i- sin(2M + k ■ 2ri„ 



sin 2rin J sin 2rin 

±6{K^{TU - 2r]n,Vn)} . ,^ \\ . I I — . (^ . 

sm(2M - 2 ■ 2r7„) ^^-^ sm(2?7„) 



i=l j=l 



^Note that in the present section we use a sHghtly different normahzation for the transfer matrix. In 
contrast to (2.4), now T{u) = Rqi(u)Ro2{u) ■ ■ ■ Ron{u). 
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Combining the fusion hierarchy with the truncation identity for a particular value n, 
yields an (n + l)-order functional equation for the transfer matrix of the model. Unfortu- 
nately, it does not seem plausible to write the functional relations in a closed form. Starting 
from low values of n, one observes that at each level new terms emerge and their number is 
given by the Fibonacci numbers. In particular, at a given value n, there are 

F(n + 1) + F(n- 1) + 1 , n = 2,3,--- 

terms in total. Nevertheless, after carefully examining the structure of the functional rela- 
tions, it is seen that the following simple schematic structure appears 

i{u) t{u + 2r])--- t{u + (2n - 2)r]) + Ui + ^l_i{u) ■ Oq 



n -1 n -I 

2" -^ 2 -^ 



lid 



{u} + JJ A(m + Akr]) + 0;_i(m) W K{u + {Ak - 2)r]) . (4.4) 

fc=0 k=l 

First, we should note that the last two terms in the RHS are present only for even values of 
n. Next, the symbol O stands for the following sequence: 



n-2 



Oi:= > A(M + 2go^)r,„(M) 



<?o=0 



ra— 4 (Ji 

+ ^K{u + (2gi + 4)r7) ^ K{u + 2^0^]) % 

<ji=0 go=0 



qmi \^) 



n-6 (?2 gi 

5^ A(m + (2g2 + 8)r]) ^ A{u + (2gi + A)r]) J] A(« + 2gor/) T,, 

92=0 gi=0 90=0 



l30'?l'?2l^J 



+ ■■■, (4.5) 



with 



n-l 



T,o9i-9.(«) = n ^~(" + 2"^^)' Q = {(li + ^h q^ + 2^ + 1}, z = 1, 2 



) ■") 



m=0 

A similar structure holds for Oq as well, although some lower/upper limits of the sums 
and the product are different. Nevertheless, the structure above is reminiscent of a path- 
ordered exponential, and appears as some deformed discrete version of the latter one. It 
would be interesting to see if this observation possesses some physical meaning, i.e., if the 
corresponding operator in the path-ordered exponential plays some physical role here. 

Despite the fact that the functional relations cannot be written in a closed form, they 
can be cast into a vanishing determinant representation, which afterwards provides the TQ 
equations. The restriction m ^ Q in (4.5) reflects exactly this property. 
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4.1 The determinant representation 



The vanishing determinant representation has the same structure with the corresponding 

one of the XXZ model [10]. In particular, the functional relations (4.4) take the following 

form: 

/ Ao -h^i ■■■ -ho\ 

-hi Ai -ho ■■■ 

-h2 h -hi ■■■ 



det 











(4.6) 



/ 



■ ■ ■ -hn-1 An-1 -hn-2 

\-hn-l ■■■ -hn K 

where A^ = K{u + 2kri) is an eigenvalue oit{u + 2kr]), hk = h{u + 2kri) and hk = h{u + 2kr]). 
The functions hk, hk are determined by demanding relations (4.4) and (4.6) to be identical. 
One concludes then that 

h{u + 21]) h{u - 21]) = -A{u) , 



h{u) h{u + 2(n - l)i]) = -^Hu) 



(4.7) 



n— 1 



lid/ 



Y[ h{u + 2ki]) + Y[h{u + 2ki]) = 4)i';[u) 

fc=0 fc=-l 

Inspired by the structure of the /i-function in the XXZ case [10], as well as the functions 
found in [18], we consider the following expressions 



h{u) = h{—u — 4ri), 



h{u) 



2N 



fsin{u + 21]) \ sin(2-u + Arj) 



9~{u)g+{u), 



{AA 



\ sin2r7 J sin(2-u + 27])' 
where the functions g±{u) contain the boundary information. Substituting the above ex- 
pressions into (4.7), the following condition arises regarding the functions g±{u): 

g_{u + 2i])g^{-u - 2i])g+{u + 2i])g+{-u - 2i]) = det K'^iu) det K~{u + 2i]) . (4.9) 

Assuming that the functions factorize the determinants of the reflection matrices as 

(7_ {u + 2ri)g^{—u — 21]) = det K~ {u + 21]) 

(4.10) 
g+{u + 2i])g^{—u — 2i]) = — det K^{u) , 

and recalling the explicit expressions for the determinants of the boundary matrices, the 
relations above are then equivalent to 

g-{u)g^(—u) = det K^{u) = — (cj^)^ sin(-u + %jj_) sin(-u — ■?/'_) 
g+{u)g^(—u) = det K^{u — 2i]) = —(cu^)'^ sm{u + ■?/'+) sin(-u — 4'+) ■, 



(4.11) 
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pinpointing to the natural solutions 

g-{u) = u' sm{u + ip_) , g+{u) = u^ sm{u + ip+) . (4.12) 

Using the expressions (4.12) then, all three relations in (4.7) are automatically satisfied. 
It should be stressed out that the construction so far does not involve the boundary odd 
Grassmann numbers at all. In fact, the structure we have found here is identical with the 
corresponding structure of the purely diagonal XXZ case [10], since the odd Grassmann 
parameters drop out from the determinants due to their nature. 

The vanishing determinant guarantees the existence of a null vector, {Qo, Qi,- ■ ■ Qn), so 
that the following relations hold 

^oQo-h-iQi-hoQn = 0, 
-hkQk-i + ^kQk-hk-iQk+i = 0, A; = l,---,n-l, (4.13) 

— hn-lQo — hnQn-1 + ^riQn = . 

For the present choice of the quasi-classical parameter, i.e. t] = rjn = ^^, these equations 
can be recast as TQ equations for a periodic function Q{u) by identifying Qk = Q{u + 2rik): 

Mu) = M«)%|^ + H-u - 2.,)5|±^ . (4.14) 

Q{u) Q[u) 

Note that the /i-functions appearing in the relation above essentially factorize the quantum 
determinant as h{u)h{—u) = (^^{2u)A{u — 21]). The derived TQ equation is actually the 
same with the one obtained directly from the fusion hierarchy, since h{u) may be identified 
with ha{u) defined in Eq. (3.8). One may proceed then with the deformation of the h- 
function, as well as that of the Q-functions to incorporate nondiagonal contributions. We 
have also checked explicitly that the deformed /i-functions (3.11), (3.18) satisfy the two 
additional constraints arising during the truncation of the hierarchy, that is Eqs. (4.7). 
Hence the two approaches converge to give the same results. 

5 The eigenstates of the model 

In the previous sections we have obtained the eigenvalues of the transfer matrix from the 
the fusion hierarchy and the resulting TQ equation. Similar as in the case of diagonal 
boundary matrices the eigenvalues can be associated to sectors labelled by the integer M of 
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parameters appearing in the solution. This is remarkable since the related U{1) symmetry 
is broken when the nondiagonal boundary conditions are applied. 

As is common to functional Bethe Ansatz approaches our solution of the spectral problem 
thus far does not provide information regarding the eigenstates of the model. However, 
exploiting the existence of the odd Grassmann numbers and their nilpotency, it is possible 
to exactly compute the M = state of the model for an arbitrary number N of chain sites, 
solely by using the derived expressions of the eigenvalues. To this end we choose /c-particle 
states. A; = 0, . . . , A^, 

\iii2---ik) =cl^cl---cljn) , l<ii<i2< ... <ik<N (5.1) 

as basis of the Hilbert space of the system, \Q) is the Fock vacuum of the system (and the 
M = eigenvector of the diagonal problem) . 

The key observation used for the construction of the unique M = eigenvector of the 
model with nondiagonal boundary conditions is that for generic chain length A^, only the 
Fock vacuum and basis states containing one or two particles contribute with a non-vanishing 
amplitude, all other sectors decouple. This allows to express the M = eigenvector of the 
nondiagonal model as 

N N N 

i=l j=l i<j 

The eigenvalue problem for the complete Hamiltonian reads as 

^1^) = (Adiag + ^ Anondiag)!^) • (5.3) 

Sphtting the terms with respect to the order of /3^ we obtain 

0{P) : (//bulk + //diag) /3" Yl ^^\'^ + ^nondiagl^) = Adiag /?" J^ b^\^) 

i i 

o{p^) : (i/buik + i/diag)/3+r $^5.,Kj) + //no„diagr 5^&rK) = 

i<j i 

= Q Anondiagl^) + Adiag P^ P' JZ ^«il^-^') " (^•^) 

Regarding the terms linear in (5^ first, after the relevant computations and splitting the 
resulting relations with respect to the appropriate excited states, one ends up with the 
following six relations that determine the coefficients hf 

/3+|l) : -t6+ + 6+[\/(Ar-2)-|cot^_+Ar+-Adiag] =0 

13 



rii): 
riiv): 

P^\i) : 



-t 6+ _i - CSC ^+ + fe^ [y (iV - 2) + ^ cot ^_-Af_- Adiag] = 

-t 62 - CSC ?/'_ + 61 [\/(iV - 2) - I cot tp_+Af+- Adiag] = 



-t &JV-1 + ^7V [^(^ - 2) + ^ cot ^_ - A/1 
-t (6t 1 + &ti) + ^^* [ViN - 3) + i cot ^_ 



— ^diagj — 

+ Ar+ - Adiag] =0, (5.5) 



where the recursion relations are vahd for 2 < i < N — 1 and can be solved analytically, 
giving 



bf 



1 



t 



At^ 



ct 



Cn 



t 



At^ 



^2 



(5.6) 



with Co = V{N — 3) + ^ cot ■ip_ +Af+ — Adiag and Cf 2 constants to be determined. Proceeding 
to the quadratic terms, the ones proportional to /3+a" and /3"a;+ first give the constraints 



CSC(^_) b^ = Anondiag = CSc(V^+) b^ . 



(5.7) 



These two constraints, combined with the set of relations (5.5) are sufficient to completely de- 
termine the constants €^2 and therefore all coefficients bf. We conclude that the coefficients 
bf are given by 



sin(V^_ + {£- l)2ri) 
sin((A^- l)2?7 + ^_ + V^+) 

sin(V^+ + {N - i)2r]) 



{5.i 



' sm{{N -l)2r] + ilj_+ilj+)' 

Concerning the rest of the quadratic terms, after some algebra and splitting the terms with 
respect to various states, we conclude that the coefficients Bke satisfy the following relations 



{E++Af_)BiN + t{BiN-i + B2n) + b+N csct/'- + bl csc^'^ 
(S+ - M+) B12 + 1 5i3 + 6+ CSC V- = 

(S2 +Af^)BN-lN + tBN-2N + &Ar_i CSC 7/^+ = 

(S4 + J\f_)BeN + t{Be,^iN + 5^+iAf + ^^jv-i) + bj esc V^+ 
(S+ - M+)Bu + t{Bu-i + Bu+i + B2i) + 6; csc^_ = 0, 
(S3 — M+jBu+i + t{Be^u^i + Bu+2) = 0, 
(Sg — J\f+)Bki + t{Bk-u + -Bfc+if + -Bfc^-i + Bki+i) = 0, 

where for the sake of readability we have defined 

= Adiag±|cot^_-y(iV-g). 







0, 1 <£< A^-1 

2<£<N 

1 <e<N -1 

l<k<N -2, 



>k 



^Q 



(5.9) 



In principle, having acquired the exact expressions for bf, the above set of relations provides 
the expressions of Bke as well. Since it is hard to solve these relations analytically, one may 
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resort to a numerical analysis for a given number of chain sites. However, solving these 
equations for small numbers of chain length, we were able to observe the emerging pattern 
which governs the coefficients B^e- In short, we have found that they are eventually given 
by the very simple expressions 

_ sin((iV-l)2r/ + ^.+^^) ,-,,-,+ x .. ^.^ 

We have explicitly confirmed that the expressions (5.10) satisfy all the generic constraints 
and relations derived above. In conclusion, the M = eigenvector of the model is completely 
determined for an arbitrary number of chain sites. 

The decoupling of higher/lower sectors, due to the nilpotency of the Grassmannian pa- 
rameters, further constraints the expressions of the eigenvectors which correspond to excited 
states. For generic values of M then, we propose that the corresponding eigenvector will be 
given by the following schematic expression 

|^)m = (ci + C2a+/3- + C3 /?+«-) |M)+C4/3+|M + l)+C5r|M+l) 

+ C6l3+l3-\M + 2) + C7a+\M -1) + C8a-\M -1) + C9a+a-\M ~2) , (5.11) 

where with \M) we denote the states with M particles present, or equivalently the states 
with M spins down in the spin picture, i.e. M excitations from the ground state. The 
number of states with the same M is given by the binomial coefficient 

MJ ~ M\{N-M)\ ' ^^'^^^ 

so that the states \M) in (5.11) are to be understood as collections of states spanning the 
degeneracy space for a particular M . In the same spirit, the coefficients q G C appearing 
also in (5.11) are to be interpreted as sets of coefficients of the degenerate states. 



6 Discussion 

In the present work, we have solved the small polaron model with nondiagonal boundary 
conditions. The eigenvalues of the transfer matrix have been extracted first by direct use 
of the fusion hierarchy and then by considering the functional relations for particular values 
of the anisotropy parameter. Both approaches yield the same set of eigenvalues, confirming 
our proposed deformations needed to account for the nondiagonal nature of the model. The 
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eigenvalues depend on two sets of Bethe roots, for which coupled Bethe ansatz equations 
have been presented. 

An interesting aspect of the model's solution is that, unlike in the case of the XXZ model 
with nondiagonal boundary conditions, no restrictions emerge for the boundary parameters. 
This extra freedom, as well as the remnants of particle number conservation leading to 
sectors of the Hilbert space labelled by the integer M appear to be inherited from the 
supersymmetric nature of the model. Furthermore, supersymmetry heavily restricts the 
structure of the eigenvectors and has rendered possible to exactly compute the M = 
eigenstate of the model. A more detailed analysis should provide the complete expressions 
of all eigenvectors. 

Since supersymmetry lifted any possible constraints between the boundary parameters 
for the small polaron, it is interesting to consider other supersymmetric models with non- 
diagonal boundaries. Particularly interesting would be an extension of the analysis of the 
supersymmetric t-J model with open boundaries [29, 30] to this case. In this case nondiag- 
onal terms breaking either the U{1) charge symmetry or the SU{2) spin symmetry of the 
model can be added. While the latter will likely lead to problems similar to those encoun- 
tered in the XXZ spin model, we expect that terms breaking the charge symmetry can be 
dealt with in a similar manner as in the small polaron model. 

Another route to follow is to consider operator valued representations of the reflection 
algebras [7, 31, 32] instead of c- number solutions, and attempt to solve the model for bound- 
ary conditions breaking the bulk symmetries. The construction of suitable, generalized 
Jordan- Wigner transformations for the nondiagonal boundary terms would provide valuable 
information in this spirit. 
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